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Abstract. The problem of the invariant classification of the orthogonal co- 
ordinate webs defined in Euclidean space is solved within the framework of 
Felix Klein's Erlangen Program. The results are applied to the problem of 
, integrability of the Calogero-Moser model. 

o : 
o ■ 

<N : 1 Introduction 

>v 

In his famous Erlangen Program [TT], Felix Klein introduced a unified point of view 
^ . according to which many different branches of geometry could be integrated into a 

single system. As is well known, this standpoint stipulates that the main goal of 
| any branch of geometry can be formulated as follows: 

> ; 

"Given a manifold and a group of transformations of the manifold, to 
O ' study the manifold configurations with respect to those features that are 

not altered by the transformations of the group." (JB|, p 67) 

: 

The term "manifold of n dimensions" in this setting describes a set of n variables 
that independently take on the real values from — oo to oo ([12], p 116). 

Motivated by this idea, one can assert that Euclidean geometry of E 3 (Euclidean 
space) can be completely characterized by the invariants of the Euclidean group 
of transformations. As is well-known, this Lie group of (orientation-preserving) 
isometries, denoted here by I(E 3 ), is a semi-direct product of the corresponding 
>- ! groups of rotations and translations. 

An important aspect of Euclidean geometry is the theory of orthogonal coordi- 
nate webs that originated in works of a number of eminent mathematicians of the 
past including Stackel [27], Bocher jl], Darboux [5j and Eisenhart JHi within the 
framework of the theory of separation of variables. Its modern developments can 
be found in the review by Benenti [2j and the relevant references therein. In partic- 
ular, it has been shown that there exist exactly eleven orthogonal coordinate webs 
which afford separation of variables for the Schrodinger and Hamilton- Jacobi equa- 
tions defined in E 3 . These coordinate webs are confocal quadrics determined by the 
Killing tensors of valence two having orthogonally integrable (normal) eigenvectors 
and distinct eigenvalues. Eisenhart's results in E 3 were extended by Olevsky [5T] to 
three-dimensional spaces of non-zero constant curvature, while Kalnins, Miller and 
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others generalized them to spaces of higher dimensions (see Kalnins ^0] and the 
references therein). 

This work is a natural continuation of the project initiated in ^J] (see also 
[T5] and where isometry group invariants and covariants of valence- two Killing 
tensors are derived and used to classify orthogonal coordinate webs of the Euclidean 
and Minkowski planes. Accordingly, we approach the problem of classification of 
the eleven orthogonal webs in E 3 from the viewpoint of the invariant theory of the 
isometry group /(E 3 ). Recall that the standard approach to the study of Killing 
tensors defined in pseudo-Riemannian manifolds of constant curvature rests on the 
fact that they can be expressed in this case as sums of symmetrized tensor products 
of Killing vectors. In contrast to the conventional view, we consider the Killing 
tensors of valence two defined in E 3 to be algebraic objects or elements of the 
corresponding vector space /C 2 (E 3 ) and define the action of /(E 3 ) in this vector 
space to derive I (K 3 ) -invariants of the valence- two Killing tensors. 

In line with the postulates of the Erlangen Program, we completely solve 
the problem of classification of the eleven orthogonal webs in E 3 by em- 
ploying the I (K 3 ) -invariants of the vector space /C 2 (E 3 ) and its subspaces. 

Our solution is based on the result of theorem 15.11 which describes the space of 
all isometry group invariants of the vector space of Killing tensors of valence two 
defined in E 3 combined with a careful study of the corresponding vector space of 
Killing vectors (the Lie algebra of the isometry group of E 3 ). 

It must be emphasized that the problem of the invariant classification of the 
orthogonal coordinate webs in E 3 is significantly more complicated than the corre- 
sponding problems in two-dimensional pseudo-Riemannian spaces of constant cur- 
vature. Apart from the obvious difficulties in dealing with a vector space of a much 
higher dimension, one has to solve the problem of the normality of eigenvectors 
of valence-two Killing tensors. More specifically, the eleven orthogonal coordinate 
webs in E 3 are generated by Killing tensors of valence two with normal eigenvectors. 
On the other hand, unlike the situation in two-dimensional spaces, in E 3 not ev- 
ery Killing tensor of valence two with distinct eigenvalues has normal eigenvectors. 
Moreover, the normality condition is equivalent to a system of non-linear partial 
differential equations (PDEs), which makes it nearly impossible to verify directly. 
The problem of finding necessary and sufficient intrinsic conditions for the eigen- 
vectors of a tensor field of valence two with pointwise distinct eigenvalues to be 
orthogonally integrable has a long history. It can be traced back to Schouten [24], 
where such conditions depending on the eigenvectors were derived. Tonolo [3U] 
subsequently determined a set of eigenvector-independent necessary and sufficient 
conditions for the Ricci tensor defined in a three-dimensional space to have normal 
eigenvectors. This criterion was shown by Schouten to be applicable to arbitrary 
(Ricci or not) valence-two tensor fields defined in an arbitrary pseudo-Riemannian 
manifold. Later, Nijenhuis [20] derived an equivalent formulation of the criterion 
introduced originally by Tonolo in terms of the components of the Nijenhuis tensor 
of the tensor field in question. In view of their respective contributions, we refer to 
these remarkable formulae throughout this paper as the Tonolo-Schouten-Nijenhuis 
(TSN) conditions and employ them to verify the normality of the eigenvectors of 
valence-two Killing tensors defined in E 3 . In addition, we determine in each case the 
coordinate transformation from the given Cartesian coordinates to the correspond- 
ing coordinate system determined by the orthogonal web. 
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As an illustration of the power of the new theory, we use it to obtain a concise 
solution to the problem of integrability of the Calogero-Moser Hamiltonian system 
defined in E 3 via orthogonal separation of variables in the associated Hamilton- 
Jacobi equation. 

The paper is organized as follows. In section we give an overview of the 
invariant theory for vector spaces of Killing tensors defined on pseudo-Riemannian 
manifolds of constant curvature. This theory is then specialized in section El to 
vector spaces of valence-two Killing tensors in Euclidean space. In section HJ we 
discuss Hamilton- Jacobi theory in the context of separation of variables and use it 
to derive canonical forms for the orthogonal coordinate webs of E 3 . The fundamental 
invariants are derived in section Eland are used in section El to classify the coordinate 
webs. Methods for transforming a given Killing tensor to canonical form are treated 
in section [7| In section El we summarize the steps in our algorithm and apply it in 
section El to determine separable coordinates for the Calogero-Moser system in E 3 . 
Finally, we draw conclusions in section EH and indicate future research directions. 

The reader will no doubt realize that our classification of the coordinate webs 
and our algorithm for determining separable coordinates for natural Hamiltonians 
in E 3 is highly computational. Nevertheless, all computations are purely algebraic in 
nature, and thus are straightforward to implement in a computer algebra system. To 
complement the paper, we have written a Maple package, called the KillingTensor 
package, which performs all the steps in our algorithm. The package is available 
through the Maple Application Centre at http:/ /www.mapleapps.com. 

2 Invariant theory of Killing tensors 

In the past decade, the classical invariant theory of homogeneous polynomials has 
become an active field of research once again (see Olver [22] and the references 
therein). The theory emerged in the nineteenth century as the intrinsic study of 
vector spaces of homogeneous polynomials under the action of the general linear 
group. Two of the authors (RGM, RGS) and Dennis The have incorporated the 
basic ideas of classical invariant theory into the study of Killing tensors defined in 
pseudo-Riemannian spaces of constant curvature under the action of the isometry 
group [UlinilinillZllIH] • This synergy of the two theories grew out of the observation 
that Killing tensors of the same valence defined in a space of constant curvature 
constitute a vector space or, more precisely, a representation space of the isometry 
group of the underlying space. Putting this observation in proper perspective allows 
one to extend the basic ideas of classical invariant theory to the study of Killing 
tensors. Indeed, let (M, g) be an n-dimensional pseudo-Riemannian manifold of 
constant curvature with metric tensor g. 

Definition 2.1. A Killing tensor K of valence p defined in (M,g) is a sym- 
metric (p, 0) tensor satisfying the Killing tensor equation 



where [, ] denotes the Schouten bracket \2S§ . When p = 1, K is said to be a Killing 
vector (infinitesimal isometry) and the equation \2. 1\) reads 



[K,g]=0, 




£k9 = 0, 



(2.2) 



where £ denotes the Lie derivative operator. 
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The Schouten bracket is a real bilinear operator, which property together with 
(j2.1j) implies that the set /C P (M) of all Killing tensors of valence p defined in (M, g) 
is in fact a vector space. Its dimension d is determined by the Delong-Takeuchi- 
Thompson (DTT) formula 01211123 

Therefore the general element of /C P (M) is represented by d arbitrary parameters 
a 1 ,..., a , with respect to an appropriate basis. Alternatively, this fact can be 
verified by solving the corresponding Killing tensor equation (j2.1j) with respect to 
a fixed system of coordinates, in which case the parameters a 1 ,..., a appear as 
constants of integration in the general form of elements of /C P (M). 

Each element h of the isometry group I(M) induces, by the push forward map, a 
non-singular linear transformation p(h) of /C P (M). By theorem 3.5 of McLenaghan 
et al [TS], the map 

p : I(M) -> GL(/C P (M)) (2.4) 

defines a representation of I{M). Indeed, p is a group isomorphism. Once the form 
of the general element K of /C P (M) is available with respect to some convenient 
system of coordinates on M, the explicit form of the transformation p(h)K (written 
more succinctly as h ■ K) may be written explicitly in terms of the parameters 
a 1 , . . . , a. We shall be particularly concerned with the smooth real- valued functions 
on )C P (M) that are invariant under the group I(M). The precise definition of such 
/(M)-invariant functions of /C P (M) is as follows. 

Definition 2.2. Let (M, g) be a pseudo-Riemannian manifold of constant curvature. 
Let p ^ 1 be fixed. A smooth function F : /C P (M) — > R is said to be an I(M)- 
invariant of JC P (M) iff it satisfies the condition 

F(h-K) = F(K) (2.5) 

for K G /C P (M) and for all h G I(M). 

The main problem of any invariant theory is to describe the whole space of 
invariants of a vector space under the action of a group. To achieve this one has to 
determine a set of fundamental invariants with the property that any other invariant 
is an analytic function of the fundamental invariants (see [221 f° r more details). The 
fundamental theorem of invariants of a regular Lie group action [22] determines 
the number of fundamental invariants needed to define the whole of the space of 
I ( M) -invariant s . 

Theorem 2.1. Let G be a Lie group acting regularly on an n-dimensional manifold 
M with s-dimensional orbits. Then, in a neighbourhood N of each point p G M , 
there exist n — s functionally independent G -invariants A 1; . . . , A n _ s . Any other G- 
invariant X defined near p can be locally uniquely expressed as an analytic function 
of the fundamental invariants through I = F(Ai, . . . , A n _ s ). 

In order to determine the form of the invariants of /C P (M), we use the fact the 
the invariance of a function under an entire Lie group is equivalent to the invariance 
of the function under the infinitesimal transformations of the group given by the 
corresponding Lie algebra. The precise result is given in the following proposition 

[221- 
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Proposition 2.1. Let G be a connected Lie group of transformations acting regularly 
on a manifold M . A smooth real-valued function F : M — > R is a G-invariant iff 

v(F) = (2.6) 

for all p £ M and for every infinitesimal generator v of G. 

In our application, G is the representation p(I(M)) denned by ()2.4j) and the 
condition (j2.fi |) is equivalent to 

I7 f (F) = 0, z = l,...,r, (2.7) 

where the Ui are vector fields which form a basis of the Lie algebra of the repre- 
sentation and r = dim I(M) = \n{n + 1). By theorem 3.5 of ^H], this Lie algebra 
is isomorphic to the Lie algebra of I{M). Such a basis may be computed directly 
as the basis of the tangent space to p(I(M)) at the identity if an explicit form of 
the representation is available. According to theorem 12. II of the present paper, the 
general solution of the system of first-order PDEs ()2.7j) is an analytic function F 
of a set of fundamental /(M)-invariants. The number of fundamental invariants is 
d — s, where d is given by ()2.3|) and s is the dimension of the orbits of p(I(M)) 
acting regularly in the space /C P (M). 

To determine s and the subspaces of )C P (M) where the isometry group I(M) acts 
with orbits of the same dimension, one can use the result of the following proposition 

Proposition 2.2. Let a Lie group G act on M and let p £ M. The vector space 
S\ p = span{C/j| p | Ui £ g} spanned by all vector fields determined by the infinitesimal 
generators at p coincides with the tangent space to the orbit O p of G that passes 
through p, i.e. S\ p = T p (O p ). In particular, the dimension ofO p equals the dimension 
ofS\ p . 

We are now prepared to apply the theory presented thus far to the vector space 
/C 2 (E 3 ). 



3 Invariant theory of Killing tensors of valence 
two in Euclidean space 

We now specialize the general theory of the previous section to the vector space 
/C 2 (E 3 ) of valence- two Killing tensors in Euclidean space E 3 . Recall the following 
well-known result in j2H| from invariant theory. 

Theorem 3.1. The orbits of a compact linear group acting in a real vector space 
are separated by the fundamental (polynomial) invariants. 

We first note that in our case the group is non-compact and so in order to 
distinguish between the orbits of /(E 3 ) acting in the vector space /C 2 (E 3 ) we need to 
employ a more elaborate analysis than a mere computation of a set of fundamental 
invariants. 

It is well-known that in E 3 , as in all manifolds of constant curvature, any Killing 
tensor is expressible as a sum of symmetrized products of Killing vectors. The six 
Killing vectors in E 3 may be written in Cartesian coordinates x % viz 

X t = A r. = e k jiX j X k , (3.1) 
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for i = 1,2,3, where is the Levi-Civita permutation tensor 4 . We also note the 
commutation relations 

[Xi,Xj] = 0, [Xi,Rj] = e k ijX k , [Ri,Rj] = e k ijR k - (3.2) 

Thus the general Killing tensor in /C 2 (E 3 ) may be expressed as 

K = A ij Xi X 3 + 2B ij Xi © Rj + C lj Ri © Rj, (3.3) 

where the coefficients A % \ B % i and C l i are constant and satisfy the symmetry prop- 
erties 

A ij = A m^ C ij = C (ij)_ (34) 

It follows from (J3.1)) and ()3.3j) that the components of the general Killing tensor in 
/C 2 (E 3 ) with respect to the natural basis are given by 

K ij = A ij + 2e^ £k B^ k x e + e l mk e\,C M x m x n . (3.5) 

For future reference, we give explicitly the six independent components of K l K 
Noting the symmetries (j3.4j) . it proves convenient to set (following [3]) 

&11 &12 &i3 N 

&21 &22 &23 I , C l] = I 73 c 2 7i I (3.6) 





( ai 


«3 


«2^ 




A ij = | 


a 3 


«2 








V«2 


O: i 


a 3 y 







fa 


73 


72 \ 


, c* = | 


73 


c 2 


7i 




V72 


7i 





. &31 &32 ^33, 

and x l = (x, y, z). From (J3.5J) we obtain 

K 11 = ai- 2b 12 z + 2b 13 y + c 2 z 2 + c 3 y 2 - 2^yz, 

K 22 = a 2 - 2b 23 x + 2b 2 iz + c 3 x 2 + c x z 2 - 2>y 2 zx, 

K 33 = a 3 - 2b 31 y + 2b 32 x + c x y 2 + c 2 x 2 - 2^ 3 xy, 

K 23 = ai + b 31 z - b 21 y + (b 22 - b 33 )x + (73Z + j 2 y - "f\x)x - c x yz, 

K 31 =a 2 + b 12 x - b 32 z + (b 33 - b n )y + (7^ + 732; - j 2 y)y - c 2 zx, 

K 12 = a 3 + b 23 y - b 13 x + (b n - 622)^ + (72Z/ + Jix - l3^)z ~ c 3^y. 



(3.7) 



According to the DTT formula ()2.3|) . the dimension of /C 2 (E 3 ) is twenty which 
appears to disagree with ()3.6|) which lists twenty-one parameters. To reconcile this, 
we observe from ()3.7|) that only the differences of the diagonal coefficients bn, b 22 
and 633 are involved. Defining 

Pi = b 22 - 633, (3 2 = 633 - 6n, (3 3 = 6n - 622, (3.8) 

yields the constraint Pi + /3 2 + fl 3 = 0, thereby showing that there are twenty 
independent parameters. In many of the computations which follow, it turns out to 
be more convenient to use the three bn parameters instead of two of the three 
With this in mind, we (commit an abuse of notation and) let /C 2 (E 3 ) be the space 
spanned by the twenty-one parameters 

ai, a 2 , a 3 , ai, "2, a> 3 , bn, b 22 , b 33 , b 23 , b 31 , b 12 , b 32 , b 13 , b 21 , c x , c 2 , c 3 , 71, 72, 73. (3.9) 



4 We are using the summation convention throughout and lowering and raising indices with the 
Euclidean metric gtj — diag(l, 1, 1) and its inverse g lJ . 
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We shall also refer to ()3.9j) as the Killing tensor parameters. 

We now consider the transformation rules for the Killing tensor parameters. The 
transformation from one set of Cartesian coordinates x % to another set x % is given 
by 

x i = X/x j + 5\ (3.10) 

where A/ G 50(3) and 5 l G IR 3 . It is straightforward to show that the Killing 
vectors (|3.1|) transform according to 



X, 



X 3 iXj, 11, = X 3 iRj + fi'jX r 



where 



(3.11) 
(3.12) 



The Killing vector transformation rules ()3.11|) in conjunction with ()3.3|) lead to 

A ke X i k X j i + 2B ke X (i k ^ j) e + C fe V 



B ij 



B M X\X\ + C Kl Xhv\, 

k' 



iki • 



(3.13) 



These equations give the explicit form of the representation of /(E 3 ) on /C 2 (E 3 ) with 
respect to a Cartesian coordinate system on E 3 . 

Equipped with these transformation rules, we can now derive the infinitesimal 
generators of J(E 3 ) in the representation defined by ()3.13|) . Let U m , m = 1,2,3, 
denote the generators associated to the Killing vectors X m . Noting that such Killing 
vectors generate translations about the x m -axis, we set Xj 1 = 5/ in (13. 13 J) and 
differentiate the resulting equations with respect to 5 m to obtain 



d5 r 



2f mkl^ i 



dB l 'J 



s i =o 



d5 r 



£ mkC^ 



dC ij 



8 i =0 



dS r ' 



0. 



^=0 



The corresponding differential operators are therefore 

U, 



dB'J k ' 



(3.14) 



for % = 1,2, 3, where the range of summation over the derivative operators is under- 
stood to be over only those parameters listed in ()3.9|) . Next, let V m , m = 1,2,3, 
denote the generators associated to the Killing vectors R m) the generators of rota- 
tions about the x m -axis. For an infinitesimal rotation about the x 3 -axis by an angle 
9 3 , the rotation A/ G 50(3) is given by 



'cos# 3 -sin# 3 X 
A/ = [ sin# 3 cos# 3 
1 



dA/ 



d6 3 







-l 






1 







3 =O 










More generally, for an infinitesimal rotation about the x m -axis, 



dA, 
dO n 



6 m =0 
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It thus follows from ()3.13j) that 



+ {e j u C lk + e k H a") 



OA* g '8B» (3 . 15) 



dCi 



jk ' 



for i = 1,2,3. As required by the general theory, the generators ()3.14j) and ()3.15j) 
satisfy the same commutation relations as the Killing vectors (13. lj) . namely 

[Ui, UA = 0, [Ui, V 3 ] = e k tJ U k , [V u VA = e k l3 V k . 

For computational purposes, we shall require explicit expressions for the generators. 
It follows from (J3~T^) and (jHUHD that 

j=l j=i J 

for i = 1, 2, 3, where a 3 , j — 1, . . . , 21, are the twenty-one Killing tensor parameters 
ordered by ()3.9j) and 



( 


-26 23 


26 32 


622 — 633 


612 




-613 





-71 


2&13 





-2631 


-621 


633 - 611 


623 


72 





-2&12 


26 2 i 





&31 


-632 




hi - b 2 2 


-73 


73 





2ai 


-2ai 


a 3 - a 2 











623 + 6 32 


-2a 2 





2a 2 


a 3 






-Oi\ -631 - &13 





V 2a 3 


-2a 3 





— «2 


a 1 




a 2 — ai 612 + 621 


-bi2 ~ 621 


7i 


-C3 


73 







c 2 




-72 


-72 







-ci 


7i 




-73 c 3 










72 







-C2 




-71 


ci 


— &23 — &32 


^33 - 


622 


-621 


6l3 


&3S 


— 6 22 — 612 


631 


&31 + &13 


&21 


611 - 633 


-632 




612 611 - 


633 


"^23 







632 622 — 611 




-63I 6 2 3 ^22 


"fell 



















\ 



























































271 


-271 c 3 


- c 2 


-73 72 








"2 7 2 





272 


73 ci 




C3 -71 








273 


-273 





-72 


71 


C2 - Ci / 


ij 





Finally, we observe that the coefficient matrix Qj has rank six almost everywhere, 
and so, in view of theorem 12 .![ we expect fifteen fundamental 7(E 3 )-invariants. The 
computation and presentation of these invariants are treated in section 

4 Hamilton-Jacobi theory and orthogonal coordi- 
nate webs 

Consider a Hamiltonian system defined on (M, g) by a natural Hamiltonian function 
of the form 

H = y ij (x)p iPj + V(x), i,j = 1, . . .n, (4.1) 
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with respect to the canonical Poisson bi-vector P = YH=i fa* A ^? given in terms 
of the position-momenta coordinates (x,p) = (x l ,pi), i — 1, . . . ,n on the cotangent 
bundle T*(M). As is well-known, in many cases the Hamiltonian system defined 
by (j4.1J) can be integrated by quadratures by finding a complete integral W of the 
corresponding Hamilton- Jacobi equation which is a first-order PDE given by 

-g l3 {x)—— + V(x) = E, Pj = -—. 4.2 

2 v ' dx l dx? v ' dx l y ' 

The geometrical meaning of equation (J4.2)) and its complete integral W is well-under- 
stood (see, for example, PQ). Thus, if the function F = ^g^W^Wj + V — E = 
is regular on the cotangent bundle T*(M), then equation (|4.2j) defines a hyper- 
surface in T*(M). Furthermore, W is a complete integral of ()4.2j) iff the La- 
grangian submanifold S C T*(M) determined by the equations Pi = lies on 
the hypersurface defined by the Hamilton- Jacobi equation (|4.2|) . Solving ()4.2j) is 
normally based on finding a canonical transformation to separable coordinates: 
(x,p) — > (u,v) with respect to which the equation can be solved under the addi- 
tive separation ansatz W(u; c) = Y^i=i Wi{u l ] c) and the non- degeneracy condition 
det(d 2 W/du l dcj) nxn ^ 0, where c = (ci, . . . , c n ) is a constant vector. Orthogonal 
separation of variables occurs in the case when the transformations to separable co- 
ordinates are point-transformations and the metric tensor g is diagonal with respect 
to the coordinates of separation (u, v). A useful criterion for orthogonal separability 
is given by Benenti 

Theorem 4.1. The Hamiltonian system defined by \4-l\ l is orthogonally separable 
if and only if there exists a valence-two Killing tensor K with (i) pointwise simple 
and real eigenvalues, (ii) orthogonally integrable (normal) eigenvectors and (Hi) such 
that 

d(-RTdV)=0. (4.3) 



A Killing tensor satisfying conditions (i) and (ii) of theorem 14. II is called a char- 
acteristic Killing tensor ( CKT). 

Let us elaborate on conditions (i) and (ii) in theorem 14.11 On two-dimensional 
Riemannian manifolds of constant curvature these conditions are trivial, since every 
eigenvector £ of K is normal and K has repeated eigenvalues iff it is a multiple of the 
metric g. In three-dimensions, E 3 in particular, the situation is far more complicated. 
Computing eigenvectors of a symmetric 3x3 tensor is tedious and becomes virtually 
intractable if one considers Killing tensors with arbitrary parameters. Instead, we 
employ the Tonolo-Schouten-Nijenhuis (TSN) conditions, as introduced in section ^ 
which are both necessary and sufficient for a given symmetric (Killing) tensor field 
to have integrable eigenvectors. These conditions read 

N'wgut = 0, (4.4a) 
N% k K i]t = 0, (4.4b) 
N\ jk K t]m K m e = (4.4c) 

where N % j k are the components of the Nijenhuis tensor of K l i given by 

N^ k = K\K\ m + K\K\ U . (4.5) 
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We remark that the TSN conditions ()4.4aj) - (j4.4cj) yield 10 quadratic, 35 cubic and 84 
quartic equations, respectively, in the Killing tensor parameters. It is thus straight- 
forward to verify, using ()4.4|) . if a given Killing tensor satisfies condition (ii) of 
theorem 14.11 However, we have been unable to solve the conditions ()4.4|) directly to 
obtain the most general Killing tensor admitting orthogonally integrable eigenvec- 
tors 5 . Condition (i), namely the distinct eigenvalues condition, like condition (ii), 
can also be verified directly for a given Killing tensor: we simply compute the dis- 
criminant of the characteristic polynomial of K l i and verify that it does not vanish 
identically. 

Although a general solution of the TSN conditions ()4.4|) appears intractable, we 
may instead employ Eisenhart's method [B] to derive all Killing tensors with normal 
eigenvectors up to equivalence. In particular, each representative Killing tensor 
characterizes separability of the Hamilton- Jacobi equation ()4.2|) in one of the eleven 
(orthogonally) separable coordinate systems in E 3 . The Eisenhart method can be 
described as follows. Consider the Euclidean metric 

ds 2 = guidu 1 ) 2 + g 2 2(du 2 ) 2 + g 33 (du 3 ) 2 , 

with respect to separable coordinates u % . The method yields three canonical Killing 
tensors given by 

Kij = diag(Ai g u , A 2 g 2 2, A 3 #33), 
where the Aj satisfy the linear system of PDEs 

d\- d 

—± = (A,, - A,) — ]ng u , (no sum) (4.6) 

for i,j = 1,2,3 (see [Hj, equations (1.8)). Trivially, any multiple of the metric g 
satisfies (J4.6|) . It is straightforward to solve ()4.6|) for each of the eleven separable 
coordinate systems in E 3 to obtain the two additional (non-trivial) canonical Killing 
tensors which we shall label K\ and K 2 - 

We now summarize the results of this calculation. For each of the eleven separa- 
ble coordinate systems, we give the corresponding coordinate transformation, ranges 
of the separable coordinates, the metric and the components of the canonical Killing 
tensors K\ and K 2 - 



Cartesian: 
(x,y,z) 
(4.7) 



x = x,y = y,z = z 
— oo < x, y, z < oo 

„2 i J„,2 i A ^ 



ds 2 



dx + dy + dz A 
diag(0,l,0) 
diag(0,0,l) 



Circular cylindrical: 

(r,9,z) 
(4.8) 



x = r cos 6, y = r sin 6, z = z 

r ^ 0, ^ 9 < 2tt, -oo < z < oo 

ds 2 = dr 2 + r 2 d8 2 + dz 2 

K\ j = diag(0,r 4 ,0) 

i^' = diag(0,0,l) 



Steve Czapor (private communication) has simplified the situation considerably. Using 
Grobner basis theory, he has shown that Ij4.4a(l and (|4.4bJ) imply (|4.4c|l . for any Killing tensor 
Ke/C 2 (E 3 ). 
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Parabolic cylindrical: 

(/!, V, z) 

(4.9) 



X = |(/i 2 - V 2 ), y = Z = Z 

/i ^ 0, — OO < V < oo, — oo < z < oo 

ds 2 = (n 2 + v 2 )(d/j 2 + dv 2 ) + dz 2 

K[ j =dmg(v 2 g 11 ,- f i 2 g 22 ,0) 

J#'=diag(0,0,l) 



Ellipt ic-hyperbolic : 

(r),ijj,z) 
(4.10) 



Spherical: 

(r, 9, 0) 
(4.11) 



x = acosh^cos-?/', y = a sinh 77 sin t/>, z = z 
77 ^ 0, ^ ip < 2tt, — oo < z < oo, a > 
ds 2 = a 2 (cosh 2 77 — cos 2 7/>)(d77 2 + dip 2 ) + d^ 2 
K\ 3 = diag(a 2 cos 2 ip gn, a 2 cosh 2 rj g 22 , 0) 
#?=diag(0,0,l) 



x = rsva.6 cos 0, y = r sin 6* sin <f>, z = r cos 6> 
r ^ 0, 0^^<tt, O^0<2tt 
ds 2 = dr 2 + r 2 d9 2 + r 2 sin 2 d0 2 
K[ j = diag(O,r 4 ,r 4 sin 2 0) 
K l 2 j = diag(O,O,r 4 sin 4 0) 



Prolate 
spheroidal: 

(4.12) 



x = a sinh r] sin 9 cos ip, y = a sinh 77 sin # sin ip, z = a cosh 77 cos 6? 
77 ^ 0, 0^9 < 7r, ^ t/> < 2vr, a>0 
ds 2 = a 2 (sinh 2 77 + sin 2 9)(dr] 2 + d9 2 ) + a 2 sinh 2 77 sin 2 d-0 2 
K\ J = diag( — a 2 sin 2 9 gn, a 2 sinh 2 7? g 22 , a 2 (sinh 2 77 — sin 2 9)g^ 
Ki j = diag(0, 0, a 2 sinh 2 77 sin 2 9 g 33 ) 



Oblate 
spheroidal: 

(4.13) 



x = a cosh 77 sin cos ip, y = a cosh 77 sin 9 sin ip, z = a sinh 77 cos 9 
77 > 0, ^ < vr, ^ t/> < 2tt, a > 
ds 2 = a 2 (cosh 2 77 - sin 2 6>)(d77 2 + d6> 2 ) + a 2 cosh 2 77 sin 2 9 dip 2 
K\ J = diag(a 2 sin 2 9 gn, o 2 cosh 2 77 g 22 , a 2 (cosh 2 77 + sin 2 9)g 33 ) 
K % 2 = diag(0, 0, a 2 cosh 2 77 sin 2 9 (733) 



Parabolic: 

(//, u, ip) 
(4.14) 



x = jiucosip, y = fiusinip, z = |(/i 2 — v 2 ) 
fi^O, i/>0, ^ ip < 2tt 
ds 2 = (/I 2 + u 2 )(dfi 2 + du 2 ) + /iV dij 2 
K{ 3 = diag( - z/ 2 #n,/i 2 #22, (/u 2 - ^ 2 )#3s) 
^ = diag(0,0,/iV(733) 
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Conical: 

(r,M) 
(4.15) 



r > 



6c J ' " 6 2 (c 2 - 6 2 ) 

6 2 < # 2 < c 2 , < A 2 < 6 2 , 



H <? = At 

1 (# 2 -6 2 )(c 2 -# 2 ) 

^ = diag(0,r 2 A 2 te ,r 2 ^33) 

#J = diag(0, r 2 g 22 , r 2 g 33 ) 



d9 2 



b 2 (c 2 ~b 2 ) 



r 2 (6 2 - A 2 ) 
(6 2 -A 2 )(c 2 -A 2 ) 



dA 2 



Paraboloidal: 
(fi, is, A) 
(4.16) 



2 _ 4(/i-6)(6-^)(6-A) 2 
b — c 

z = fi + u + \ — b — c 
Q<u<c<\<b<n<oo 



4(/i - c)(c- z/)(A - c) 



6-c 



ds 2 



(fJ--v)(f*- A) 
(li-b)U-c) 



d/x 2 + 



dA 2 



(/i- z/)(A - z/) 
(6-i/)(c-i/) 



dz/ 2 



(6-A)(A-c) 
= diag(2(z/ + A)<7n, 2(A + /i)# 22 , 2(/i + z/)# 33 ) 
j i , -4A/i£f 2 2, -k\ivg 33 ) 



i^ 2 J = diag(-4i/A^f 



Ellipsoidal: 

(77,0, A) 
(4.17) 



2 . (o-iy)(a-g)(o-A) 2 _ (6 - r,)(b - 9){b - A) 
x , , „ , \ ,1/ 



(a — 6) (a — c) 
z2 (c - V )(c - 9)(c - \) 

(c — a)(c — 6) 

a>7]>b>8>c>\ 

(V - 9)(V - A) 



ds 



d?7 2 + 



(6-a)(6-c) 



A{a-rj){b-rj)(c-rj) A(a-9)(b-9)(c-9) 
(A - ??)(A - 9) Jx2 



dA 2 



4(a- A)(6- A)(c-A) 
isf = diag( - (9 + A)# u , -(A + 77)^22, -(77 + % 33 ) 
k K l 2 ° = diag(9Xg 11 , Xr]g 22 ,r]9g 33 ) 



As we are dealing with Hamiltonian systems defined in terms of Cartesian co- 
ordinates, the next step is to transform the components of each of the canonical 
Killing tensors K\ 3 and K l 2 3 to Cartesian coordinates. This again is a routine cal- 
culation using the transformations from Cartesian to separable coordinates listed 
in (|4. 7)1 - 1)4.1 7)1 and the appropriate tensor transformation law. For each of the 
eleven separable cases, we take a linear combination of K[ 3 , K 2 and the metric 
g % K Using (J3.7j) . we then identify the constants in the linear combination and any 
essential parameters 6 with the Killing tensor parameters (J3.9)) appearing in (J3.7I) . 
We note that if the separable case under consideration has n essential parameters, 
then one can generally choose n + 3 of the Killing tensor parameters in the iden- 
tification. However, in the paraboloidal and ellipsoidal cases, it is convenient to 



6 These refer to any parameters appearing in the canonical Killing tensors K 1 ^ and iC£ ■ 
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choose more than n + 3 parameters in the identification so that the components of 
the resulting Killing tensor are polynomials in the parameters and the Cartesian 
coordinates. Consequently, this leads to algebraic constraints in the Killing tensor 
parameters. These constraints not only ensure that the resulting Killing tensor has 
normal eigenvectors, but also guarantees that one can always (uniquely) recover the 
original constants in the linear combination and all essential parameters from the 
identified Killing tensor parameters. Each of the CKTs constructed in this manner 
uniquely defines one of the eleven possible orthogonal coordinate webs in E 3 . We 
now present the results of this procedure. For each of the eleven separable coordi- 
nate systems, we give the components of the corresponding CKT with respect to 
Cartesian coordinates and any restrictions on the Killing tensor parameters. 

1. Cartesian web 




(4.18) 



2. Circular cylindrical web 




-c 3 xy 
ai + c 3 x 2 





(4.19) 



3. Parabolic cylindrical web 




ai - 2b 23 x 





(4.20) 



4. Elliptic-hyperbolic web 




(4.21) 



5. Spherical web 




-c 3 xy 
+ c 3 x 2 + c 2 z 2 

-c 2 yz 




(4.22) 



6. Prolate spheroidal web 




-c 3 xy 
+ c 3 x 2 + c 2 z 2 
-c 2 yz 




(4.23) 



7. Oblate spheroidal web 




-c 3 xy 
+ c 3 x 2 + C 2 Z 2 
~c 2 yz 




(4.24) 
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8. Parabolic web 



K i3 



ai - 2b 12 z + c 3 y 2 
-c 3 xy 
b 12 x 



-c 3 xy b 12 x\ 
a x - 2b 12 z + c 3 x 2 b 12 y 
h 2 y ai j 



(4.25) 



9. Conical web 



K 



■1 j 



'a x + c 2 z 2 + c 3 y 2 
-c 3 xy 
—c 2 zx 



-c 3 xy 

di + C3X 2 + C\Z 2 

-ctyz 



—c 2 zx 
-ciyz 
ai + ciy 2 + c 2 x 2 



(4.26) 



10. Paraboloidal web 



la x - 2b 12 z + c 3 y 2 



\ 



-c 3 xy 
b\ 2 x 



-c 3 xy b 12 x \ 

a 2 + 2b 21 z + c 3 x 2 -b 21 y 
-b 2 \y a 3 j 



bi 2 [bi 2 b 21 + c 3 (a 2 - a 3 )} + b 21 [b 12 b 21 + c 3 (a x - a 3 )] = 



(4.27) 



11. Ellipsoidal web 



'a x + c 2 z 2 + c 3 y 2 -c 3 xy -c 2 zx 

K % i = I —c 3 xy a 2 + c 3 x 2 + c\z 2 

—c 2 zx —ciyz a 3 + c\y 

(01 — a 2 )c\c 2 + (a 2 — a 3 )c 2 c 3 + (a 3 — a\)c 3 ci = 



ciyz 
2 1 c 2 x 2 



(4.28) 



We remark that the eleven Killing tensors ()4.18|) - (j4.28|) represent all possible 
Killing tensors with normal eigenvectors, up to equivalence. We have essentially 
computed the general solution of the TSN conditions (|4.4jl using Eisenhart's method. 

As we shall see in sectional the fundamental invariants of /C 2 (E 3 ) fail to discrim- 
inate amongst the eleven coordinate webs. In anticipation, we make the following 
key observation. The eleven orthogonal coordinate webs can be divided into three 
groups according to table More precisely, we say that a CKT K 6 /C 2 (E 3 ) 
is translational (rotational) if it admits a translational (rotational) Killing vector 7 
V G /C 1 (E 3 ) satisfying CyK = 0. This definition still lacks complete precision for 
we have not defined translational and rotational Killing vectors. Certainly, one can 
give a definition of such Killing vectors in terms of integral curves. But ideally we 
would like to give a definition in terms of algebraic invariants. We will revisit this 
problem in section lfi.ll For now, it suffices to note that the canonical translational 
CKTs (I4.18j) - (j4.21j) admit the Killing vector V = X 3 , while the canonical rotational 
CKTs (|Q2 j) - (jQ5|) admit the Killing vector V = R 3 . 

We now pose the following problem: construct a subspace of /C 2 (E 3 ) consisting 
of translational Killing tensors, say with V = X 3 . To proceed we take the general 
Killing tensor from (jH.H)) and impose the condition Cx 3 K = 0. This gives a linear 



7 The circular cylindrical tensor Ij4.19|l also admits a rotational Killing vector and can therefore 
be considered as both translational and rotational. 
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Table 1: The orthogonal coordinate webs in Euclidean space. 



Translational webs Rotational webs Asymmetric webs 

Cartesian spherical conical 

circular cylindrical prolate spheroidal paraboloidal 

parabolic cylindrical oblate spheroidal ellipsoidal 

elliptic-hyperbolic parabolic 



system of equations in the Killing tensor parameters which can be readily solved to 
yield 

(ai + 2bx 3 y + c 3 y 2 a 3 - b 13 x + b 23 y - c 3 xy a 2 - fiiy\ 
oi 3 - b 13 x + b 23 y - c 3 xy a 2 - 2b 23 x + c 3 x 2 ai + (3\x . (4.29) 
a 2 - fixy ati + fixx a 3 J 

We see that the four canonical translational CKTs are all special cases of ([4.29)1 . 
However, it follows that the Killing tensor (J4.29jl does not generally have normal 
eigenvectors. Consequently, we cannot take our subspace to be those Killing tensors 
of the form ([4.29)1 . Nevertheless, using the TSN conditions (|4.4jl . it can be shown that 
a.i = a 2 = (3i = is a sufficient condition for ([4.29)) to have normal eigenvectors. 
This is not a necessary condition, for the constant Killing tensor K % i = is 
of the form ([4.29)1 and has normal eigenvectors. But, it is immediate from the 
transformation rules (J3.13jl that a CKT is Cartesian if and only if it is a constant 
Killing tensor. Moreover, it can be shown from ([4.4)1 that if a Killing tensor of the 
form ([4.29)1 has normal eigenvectors and is not Cartesian, then a.\ — a 2 — fti — 0. 
Therefore, any translational Killing tensor with orthogonally integrable eigenvectors 
which is not Cartesian has the form 

(ai + 2b 13 y + c 3 y 2 a 3 - b 13 x + b 23 y - c 3 xy \ 
"3 - 6i3£ + b 23 y - c 3 xy a 2 - 2b 23 x + c 3 x 2 . (4.30) 
a 3 ) 

We define the subspace /Cf,(E 3 ) of /C 2 (E 3 ) to be the set of all Killing tensors of the 
form ([4.30)1 . We shall refer to this subspace as the space of translational Killing 
tensors, bearing in mind that it does not include all of the Cartesian CKTs. This 
does not pose any obstacle whatsoever in our goal of constructing a classification 
scheme, as we pointed out that Cartesian CKTs are trivially the constant Killing 
tensors. Finally, we remark that the space /C|.(E 3 ) enjoys two nice features. Firstly, 
the form of the general translational Killing tensor ([4.30)1 is invariant under the 
(orientation-preserving) isometry group /(E 2 ), and secondly, the upper 2x2 block of 
([4.30)1 is the general Killing tensor on the Euclidean plane /C 2 (E 2 ) (see, for example, 
[TB]). Consequently, we can take advantage of known results in the literature for 
classifying the translational webs. 

We can perform a similar analysis for the rotational webs. It follows that any 
Killing tensor with normal eigenvectors admitting a Killing vector V = R 3 has the 
form 

at — 2b\ 2 z + c 2 z 2 + c 3 y 2 —c 3 xy b\ 2 x — c 2 xz 

K% = | -c 3 xy ai - 2b 12 z + c 3 x 2 + c 2 z 2 b\ 2 y - c 2 yz 

b\ 2 x - c 2 xz bi 2 y - c 2 yz a 3 + c 2 x 2 + c 2 y 2 / 

(4.31) 
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We define the subspace JC 2 R (E 3 ) of /C 2 (E 3 ) to be the set of all Killing tensors of 
the form (j4.31J) and shall refer to this subspace as the space of rotational Killing 
tensors. We remark that the form of the general rotational Killing tensor (j4.31|) is 
also invariant under the isometry group (i.e. the group of translations about 
the z-axis) and that all canonical rotational CKTs (|4.22jl - (|4.25j) are special cases of 
(JOU). 

In conclusion, we have defined the following four vector spaces: /C 1 (E 3 ), /Cf.(E 3 ), 
/C|,(E 3 ) and /C 2 (E 3 ). For each of these four spaces, we need to compute the funda- 
mental invariants under the action of the appropriate isometry group and classify 
the corresponding canonical forms. This is the topic of the next two sections. 

5 Fundamental invariants of Killing tensors in Eu- 
clidean space 

In this section we derive the fundamental invariants in each of the four vector spaces 
/(^(E 3 ), /C^E 3 ), /C|.(E 3 ) and /C 2 (E 3 ), under the action of their corresponding isom- 
etry group. 

5.1 The space of Killing vectors 

The most general Killing vector V G /C 1 (E 3 ) may be expressed as 

V = A i X i + C i R i , (5.1) 

where Xi and Ri, i = 1,2,3, are the Killing vectors defined in ()3.1|) and the coeffi- 
cients A 1 and C l are constant. For sake of convenience, we set 

A 1 = (01,02,03), C l = (ci,c 2 ,c 3 ), 

so that the space of Killing vector parameters is spanned by the six parameters 

oi,a 2 ,a 3 ,ci,c 2 , c 3 . (5.2) 

For future reference, we note that can be written in the form 

V — (01 - c 2 z + c 3 ?/)X 1 + (a 2 - c 3 x + c l z)X 2 + (o 3 - c x y + c 2 x)X 3 . (5.3) 

As in section OJ we can derive the transformation rules for the Killing vector param- 
eters under the action of J(E 3 ). It follows that 

A* = A i X i j + CVj, & = C'X^, (5.4) 

which lead to the following infinitesimal generators of /(E 3 ) in the representation 
defined by (jOI) : 

Ui = e j ik C k ^- } Vi = e j ki A k ^- + e j ki C k ^-, (5.5) 

for % = 1,2,3. Using the formalism of section it follows from (|5.5j) that /C X (E 3 ) 
admits two fundamental J(E 3 )-invariants, namely 

Ax = C'Q, A 2 = A'Ci. (5.6) 

In the next section, we will use these invariants to classify the elements of /C X (E 3 ) 
and define translational and rotational Killing vectors in terms of Ai and A 2 . 
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5.2 The space of translational Killing tensors 



In section |3] we denned the space of translational Killing tensors /C^(E 3 ) to the set 
of Killing tensors of the form (|4.30jl . We pointed out that the upper 2x2 block 
of (J4.30)) is the general Killing tensor in the space /C 2 (E 2 ). The fundamental J(E 2 )- 
invariants of this vector space are known (see, for example, an d hence are 

also J(E 2 )-invariants of /Cf.(E 3 ). It follows that /C^(E 3 ) admits two fundamental 
J(E 2 )-invariants given by 8 

A 1 = c 3 , A 2 = [6 1 3 2 -6 2 3 2 + C3(a 2 -ai)] 2 + 4(6 13 6 2 3-a3C3) 2 (5.7) 
(see proposition 4.1 in [To]). 



5.3 The space of rotational Killing tensors 

The set of all Killing tensors of the form ()4.31j) defines the space of rotational Killing 
tensors /C|,(E 3 ). This subspace of /C 2 (E 3 ) is mapped to itself under the action of the 
isometry group /(K), the group of translations about the z-axis. Trivially, the Lie 
algebra z(R) is generated by the Killing vector field X% and hence the corresponding 
infinitesimal generator in the parameter space is the generator U3 in (I3.1fi|) restricted 
to /C|(E 3 ) which reads 

U 3 = -2b 12 —-c 2 —. (5.8) 

(7(2 x <70 12 

Solving the PDE U^F) = by the method of characteristics, we obtain the four 
fundamental /(IR)-invariants of /C^(E 3 ), namely 

A 1 = e 2 , A 2 = b 12 2 + c 2 (a 3 - at), A 3 = a 3 , A 4 = c 3 . (5.9) 



5.4 The space of Killing tensors 

We now briefly describe how to derive a complete set of invariants for the full vector 
space /C 2 (E 3 ) of valence- two Killing tensors in Euclidean space. As mentioned at 
the end of section El /C 2 (E 3 ) admits fifteen fundamental J(E 3 )-invariants which can 
be computed by solving the system of linear PDEs 

Ui(F) = 0, Vi(F) = 0, z = l,2,3, (5.10) 

where the generators U i and Vi are given by 1)3.16)1 and F is an analytic function 
in the Killing tensor parameters. 

Computationally, arriving at the general solution of (J5.1U)) is non-trivial. In par- 
ticular, we have found that the method of characteristics becomes intractable when 
applied to (|5.10j) . However, we have successfully computed all fifteen invariants us- 
ing the method of undetermined coefficients (SJ. The simplest implementation of 
this method is to build monomial trial functions in the Killing tensor parameters up 
to a fixed degree, take a linear combination of these monomials and substitute the 
combination into ()5.10|) leading to a large (sparse) system of linear equations in the 



8 The paper ^H] only treats the space of "non-trivial" Killing tensors in /C 2 (E 2 ). The invari- 
ants H5.7f) form a complete set of fundamental J(E 2 )-invariants for this vector subspace. It can be 
shown that ZC^(E 3 ) admits two additional fundamental 7(E 2 )-invariants; we do not present them 
here as they play no role in classifying the elements of JC^(K 3 ). 
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undetermined coefficients. This approach has two obvious disadvantages. Firstly, 
such an ansatz does not take advantage of the apparent structure and symmetry 
of (|5.10p . Consequently, as one must construct all monomials up to and includ- 
ing degree five to recover all fifteen fundamental invariants, approximately 50 000 
undetermined coefficients are involved; the corresponding linear system requires al- 
most ninety hours of CPU time to solve on a modest Sun workstation! Secondly, 
by virtue of the ansatz, the computed invariants are in expanded form and occupy 
fifteen pages of output. A more effective ansatz involves constructing scalar trial 
functions which are "tensorial" 9 in the A 1 - 7 , B % i and C l K For example, trial functions 
which are cubic in the include 

i ^k i ^ ^ij^k i ^ tyj ^kii t ikm t jlrC^' ^ ^ 

Implementing the method of undetermined coefficients with this ansatz yields the 
following fifteen fundamental J(E 3 )-invariants of /C 2 (E 3 ): 

r>;. A 2 c;. a, \YH-, r A; r"r, r A r , irn„ ■ A^c, r 

B^C^Cu, A 7 = C»Cj k C ki , A 8 = C»[Bj k (B ik + 2B ki ) + A k C kl \, 
e ikm e j£n B»B kl B mn - 2{B^B^ + A^C tJ )B k k + 6B ij Aj k C ki , 
B v \B k C kj - 2B k C kl ) - iir-'P,,, + A^C l3 )C k k + AiC^C k k \ 

ryii jykt r~imn r~i p , jyij r jy /nk£/n si kir~i I ry i a /"» t jy \1 
£i£m£jhp£> & ^ ^"n + -D [DijLy U k g — Uj [U k tS it + 4G[ fc rS£\i)\ 

^A^C l3 C k ^ k Cl\ 
Ai%CjC k h + 3C jk C jk )C/ - AC' k C k e C ej ] - 6A^C tJ C M C ke 
+ QB l >{B l0 C M C M - C k [{B lk - 2B ki )C/ + AC k e B £l ]} 

i in, , ryij jyki rirari r~i p 
~r i -*Htm t jkp J - > D L/ (_y n , 

/l*jVo C \kr~i (\ I jy kri ir~i n/~< ry kr~ii\ 

A J {£>ijU k L U£ + tfj U k L>n — ZOj(ji3 fc ) L>t J 

+ A i t C jk (Bj k C/ — B k Cij) — B l i[BijB ke C k £ + 2C k B ki Bi 
+ Bj k B ik C/ — (Bj k Ba + B k Btj)C k \ 
4A [ M/' 1 C fe [ fe C/] + ^{AfC^Ctf + A k k C { fC t]i ) + A^C i3 {A M C ki 
+ 4B k [k B/^) + ACFBfAk'Bu + l^CfB^B^ 

A Ojj^Ojj Of -oO Ly k e)Ly m + ZO 0£ O m fcJ 

C /iii /~i ks~1 /~i ml i r\fiij jy k { /~i try m , o D ^/"^ m\ 

— bA J Uj Uki^e — J-^O -Dj [yk tfi^m] + ^-Dfc <^i[^m] )■ 

(5.11) 

We have hence proven the following theorem. 

Theorem 5.1. Consider the vector space /C 2 (E 3 ). Any algebraic I invariant 
X defined over /C 2 (E 3 ) in terms of the Killing tensor parameters hS. .9)) where the 
isometry group /(E 3 ) acts freely and regularly with six- dimensional orbits can be 
locally uniquely expressed as an analytic function X = F(Ai, . . . , A 15 ) ; where the 
fundamental I (W?) -invariants A i; i — 1, . . . , 15, are given by H5.ll]) . 

We close this section by addressing a minor computational issue. It is clear 
from f)3.7j) that the Killing tensor parameters b\\, 622 and 633 are not uniquely 



9 As the A i: > and B y do not transform as tensors (see H3.13(l ). there is no reason why this ansatz 
should work. 



A 6 = 
A 9 = 

A10 = 
A n = 

A12 = 



A 13 = 



A 14 = 

A15 = 
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Table 2: Invariant classification of Killing vectors in Euclidean space. 



Classification 


Invariants 






translational 


A 1 = 0, 


A 2 


= 


rotational 


Ai^O, 


A 2 


= 


helicoidal 




A 2 


7^0 



determined for a given Killing tensor K G /C 2 (E 3 ). Thus, how does one evaluate the 
invariants ()5.11j) on a given Killing tensor? Indeed, we require an invariant method 
for solving the equations ()3.8j) for the b„ in terms of the known ft. This problem 
is easily rectified upon observing that the equations (|3.8|) in conjunction with the 
condition A x = <^ fen + fe 22 + fe33 = yields a unique solution for the ba given by 

fell = |(ft - ft), &22 = f (ft - ft), fe 33 = |(ft - ft)- (5.12) 

In what follows, we shall extract the parameters fen, fe 22 and &33 from a given K G 
/C 2 (E 3 ) using (EH2T) . 

6 Invariant classification of orthogonal coordinate 
webs in Euclidean space 

In order to build a classification scheme for the orthogonal coordinate webs in Eu- 
clidean space based on the set of Killing tensors in /C 2 (E 3 ) with normal eigenvectors 
and distinct eigenvalues, one must first know how to classify elements in the vector 
spaces /C 1 (E 3 ), /C^(E 3 ) and /C^(E 3 ). The classification of elements in these spaces 
are treated in subsections 16.11 16.21 and 16.31 respectively. We then use these results 
to classify the orthogonal coordinate webs of /C 2 (E 3 ) in subsection 16.41 

6.1 Classification of Killing vectors 

We shall classify the elements of /C 1 (E 3 ) according to whether the fundamental 
invariants ()5.6|) are zero or non-zero. There are three cases to consider each of which 
gives rise to a canonical Killing vector. The classification scheme is summarized in 
table El 

1. Ai = 0, A 2 = 0. In this case, the Killing vector V is of the form 

V = A l X % , (6.1) 

where it is assumed that the A % are not all zero so that V is non-trivial. It 
follows that we can use the isometry group /(E 3 ) to transform ()6.1)) to 

V = a 3 X 3 , (6.2) 

for some S3 7^ 0. Indeed, the transformation rules (|5.4jl reduce to C % = 
and A % = A : 'X t j. Without loss of generality, we can set the components of 
the translation 5 l = 0. The components of the rotation A/ can be computed 
by setting A 3 j = (A k A) s )~ 1 / 2 Aj and then obtaining A 1 ^ and A 2 j by extending 
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Table 3: Invariant classification of translational Killing tensors in Euclidean space. 



Orthogonal coordinate web 


Invariants 








Cartesian 


A! 


= 0, 


A 2 


= 





circular cylindrical 


Ai 


7^0, 


A 2 







parabolic cylindrical 


Ai 


= 0, 


A 2 







elliptic-hyperbolic 


Ai 


7^0, 


A 2 








the vector A 3 -,- to a proper orthonormal basis in E 3 (using the Gram-Schmidt 
algorithm or QR decomposition, for example). We observe that (J6.2)) defines 
a translation, hence we say that a Killing vector V G /C 1 (E 3 ) is translational 
iff Ax = A 2 = 0. 

2. Ai^O, A 2 = 0. Using the isometry group /(E 3 ), we claim that any such 
Killing vector V can be transformed to 

V = d 3 R 3 , (6.3) 

for some c 3 ^ 0. Indeed, we can first make a translation x l = x 1 + 5 l so 
that A 1 = in the new coordinates x % . It follows from the transformation 
rules (15. 4 j) and ()3.12|) that the 5 l must satisfy the system of linear equations 
e' l jkC^5 k = A 1 , which has a solution iff A 2 = 0. In the new coordinates x l , it 
follows that V = C l Ri. By a similar argument to that used in the previous 
case, we can find a rotation A/ such the coordinate transformation x 1 = A/a;- 7 
puts V into the form ()6.3|) . We observe that (|6.3|) defines a rotation, hence 
we say that a Killing vector V G /C 1 (E 3 ) is rotational iff A x ^ and A 2 = 0. 

3. Ai 7^ 0, A 2 7^ 0. By similar arguments to those used in the previous two cases, 
it can be shown in this case that there exists a coordinate transformation (j3.10|) 
such that 

V = 5 3 X 3 + c 3 R 3 , (6.4) 

for some non-zero a 3 and c 3 . The integral curves of such a Killing vector field 
are helices and so we say that a Killing vector V G /C X (E 3 ) is helicoidal iff both 
Ai and A 2 ^ 0. We remark that we can further refine this classification 
into left and right-handed helicoidal Killing vectors according to the sign of 
A 2 . 

6.2 The space of translational Killing tensors 

A classification scheme for the vector space of translational Killing tensors /C|i(E 3 ) 
based on the fundamental 7(E 2 )-invariants ()5.7|) is provided in [T5] . For complete- 
ness, we summarize this scheme in table 01 

6.3 The space of rotational Killing tensors 

Evaluating the /(IR)-invariants (IBUl) on each of the rotational CKTs 

produces a classification scheme for the vector space /C|(E 3 ). It turns out that we 
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Table 4: Invariant classification of rotational Killing tensors in Euclidean space. 



Orthogonal coordinate web Invariants 



l^LLlcLl 1111U.1 IL-dl 


^1 




n 


An 




U 


spherical 


Ai 




o, 


A 2 







prolate spheroidal 


Ai 




o, 


A 2 


> 





oblate spheroidal 


Ai 




o, 


A 2 


< 





parabolic 


Ai 




o, 


A 2 








only need to use the invariants Ai and A 2 in (|5.9jl to obtain a classification. As 
we pointed out in section EJ we can also include the circular cylindrical web in this 
classification. Our classification scheme is detailed in table EJ 

6.4 The space of Killing tensors 

The motivation for constructing invariant classification schemes in the vector spaces 
treated in the previous three subsections is due to the fact that we have been unable 
to obtain a scheme based solely on the fifteen fundamental J(E 3 )-invariants of the 
full space /C 2 (E 3 ) presented in 1)5.11)1 . This is primarily because these invariants fail 
to discriminate amongst some of the canonical CKTs (J4.18j) - (J4.28j) . 

To begin, let K e /C 2 (E 3 ) be the given CKT for which we wish to classify. 
As we showed in section 0J if K is constant, then it necessarily characterizes a 
Cartesian web. Let us therefore assume for the remainder of this section that K is 
not constant. The classification of K involves two main steps: 

1. Determine whether K characterizes a translational, rotational or asymmetric 
web (according to the type of Killing vector it admits). 

2. Use the classification schemes in tables El and 0] if K is translational or rota- 
tional, or, the classification scheme outlined in table El if K characterizes an 
asymmetric web (i.e. K admits no Killing vector). 

To proceed with the first step, we let V be the general Killing vector from (J5.3|) 
and impose the condition 

C V K = 0. (6.5) 

Equation (Jd5j) results in a linear system of equations in the six Killing vector pa- 
rameters (J5.2j) which can be readily solved. It follows that the general solution of 
()6.5)) can be decomposed as 

V = £ 1 V 1 + ---£ n V n , 

for some n ^ 6, where £j, i — 1, . . . , n, are arbitrary non-zero constants and { Vi, . . . , 
V n } is a linearly independent set of Killing vectors. If n — 0, we conclude that 
K does not admit a Killing vector, and hence characterizes an asymmetric web. 
Otherwise, using table El we classify each of the Vi according to whether they are 
translational, rotational, or helicoidal. Therefore, if one of the Vi is translational, 
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then K characterizes a translational web, otherwise K characterizes a rotational 
web 10 . 

We have now shown how to determine if K characterizes a translational, rota- 
tional or asymmetric web. To proceed, suppose that K is a translational (rotational) 
Killing tensor and let V be its corresponding translational (rotational) Killing vec- 
tor. From the results of subsection 16.21 ()6. 3)) . we can use the isometry group /(E 3 ) to 
bring V to the canonical form a^X 3 (C3-R3). Applying the corresponding coordinate 
transformation to the Killing tensor K places it in the subspace /C^(E 3 ) (/C^(E 3 )). 
Finally, we can classify the transformed K using table |3] ©. 

Suppose now that K characterizes an asymmetric web. Using the fundamental 
J(E 3 )-invariants (|5.11jl of /C 2 (E 3 ), we shall derive a scheme for classifying K. To 
begin, we evaluate the invariants A2, A4 and A7 (see (jft.llj) ) on the three asym- 
metric CKTs (|4.26J) - (|4.28|) . It follows that for the conical and ellipsoidal tensors, 
(A 2 , A 4 , A 7 ) = (ci + c 2 + c 3 , C1 2 + c 2 2 + c 3 2 , C1 3 + c 2 3 + c 3 3 ), while for the paraboloidal 
tensor, (A 2 , A 4 , A 7 ) = (03, C3 2 , C3 3 ). This motivates defining two auxiliary invariants 

E x = A 2 2 - A 4 , ~ 2 = A 2 3 - A 7 , (6.6) 

noting that S 4 = S 2 = on the paraboloidal tensor (|4.27J1 . We claim that the van- 
ishing of and S 2 is also a sufficient condition for K to characterize a paraboloidal 
web. Indeed, it follows that Hi = H 2 = in one or more of the following three cases 

ci = c 2 = 0, c 2 = c 3 = 0, c 3 = ci = 0. 

The conical case ()4.26j) and the ellipsoidal case 1)4.28)1 cannot have c\ = c 2 = 0, since 
this condition reduces them to the spherical and elliptic-hyperbolic tensors, respec- 
tively. Moreover, by a rotation, the two other cases are also impossible for conical 
and ellipsoidal tensors. Therefore, we conclude that K characterizes a paraboloidal 
web if and only if Hi = H 2 = 0. 

Suppose now that K does not characterize a paraboloidal web for the remainder 
of this section. It is convenient to define 

S 3 = 3A 4 - A 2 2 , (6.7) 

noting that H 3 = (c\ — c 2 ) 2 + (c 2 — C3) 2 + (03 — ci) 2 on both the conical and ellip- 
soidal CKTs (j4.26j) and (|4.28jl . Indeed, if H 3 = 0, the web cannot possibly be a 
conical tensor since c\ = c 2 = c 3 reduces it to a special case of the spherical tensor. 
Therefore, if H 3 = 0, then K necessarily characterizes an ellipsoidal web. 

To distinguish between a conical and ellipsoidal tensor in the case when H 3 ^ 0, 
we define three additional auxiliary invariants given by 

H 4 = A 2 A 5 - 3A 8 - 2Ai , 

H 5 = A 2 Ai + A 4 A 5 - An, ( g) 

S 6 = A 2 [2A 2 (10A 2 A 5 + 24A 8 - 3A 10 ) - 72A U + A 12 ] 
- 48A 4 A 8 - 20A 5 A 7 + I6A15 



10 It is impossible for all of the V{ to be helicoidal. Although the circular cylindrical web is the 
only coordinate web which admits a helicoidal Killing vector, it also admits both translational and 
rotational Killing vectors. Clearly, if one of the V t is helicoidal, we can conclude immediately that 
K characterizes a circular cylindrical web. 



22 



Table 5: Invariant classification of asymmetric Killing tensors in Euclidean space. 



Orthogonal coordinate web Invariants 

paraboloidal (S 1; S 2 ) = (0, 0) 

ellipsoidal (Si, S 2 ) ^ (0, 0), S 3 = or 

(Si,S 2 )^(0,0), S 37 ^0, (S 4 ,S 5 ,S 6 )^ (0,0,0) 

conical (~i,~ 2 ) + (0,0), S 3 + 0, (S 4 , S 5 , S 6 ) = (0, 0, 0) 



It follows that the invariants ()6.8|) all evaluate to zero on the conical tensor ()4.26|) . 
We claim that S 4 = S5 = = (in conjunction with S3 7^ 0) is also a sufficient 
condition for a conical tensor. Indeed, arguing by contradiction, it follows that 

5 4 = (ai + a 2 - 2a 3 )cic 2 + (a 2 + a 3 - 2ai)c 2 c 3 + (a 3 + ai — 2a 2 )c 3 ci, 

5 5 = (cic 2 + c 2 c 3 + c 3 ci)[(ai + a 2 - 2a 3 )c 3 + (a 2 + a 3 - 2ai)ci 

+ (a 3 + ai - 2a 2 )c 2 ], 

5 6 = 12cic 2 c 3 [(2ai - a 2 - a 3 )ci + (2a 2 - a 3 - a 1 )c 2 + (2a 3 - a x - a 2 )c 3 ], 

on the ellipsoidal tensor ()4.28|) and that S4 = S5 = = in one or more of the 
following five cases 11 : 

ci = c 2 = 0, c 2 = c 3 = 0, c 3 = ci = 0, ci = c 2 = c 3 , a x = a 2 = a 3 . 

By previous arguments, the first three cases are impossible if the tensor is ellipsoidal. 
The fourth case can also be eliminated since S 3 7^ 0. Finally, the fifth case is 
impossible since it reduces to a conical tensor. Therefore, we conclude that if S 3 7^ 0, 
then the tensor is conical if and only if S4 = S5 = = 0. 

This completes the derivation of the classification scheme for the set of asym- 
metric CKTs in Euclidean space. These results are summarized in table 03 



7 Transformations to canonical form 

Once a CKT in Euclidean space has been classified using the scheme detailed in 
the previous section, the isometry group J(E 3 ) can be used to transform the Killing 
tensor into its corresponding canonical form. This step leads directly to the transfor- 
mation to separable coordinates. In this section we provide methods for determining 
the transformations to canonical form. As we shall see, the majority of the calcula- 
tions amount to elementary linear algebra. 

The procedure can be summarized as follows. Suppose that are the compo- 
nents of a CKT with respect to Cartesian coordinates x % . Under the action of /(E 3 ), 
the transformation from the original set of Cartesian coordinates x l to another set 
x l is given by ()3.10j) . i.e. 

x l = \/x j + S\ (7.1) 

Thus, we need to determine the rotation A/ G SO(3) and translation 5 l G M. 3 which 
brings K % i to its appropriate canonical form given by one of the eleven cases 



1 This calculation is facilitated by use of the Maple Grobner basis package. 
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(|4.18jl - (|4.28jl (in the coordinates x l ). Moreover, any essential parameters appear- 
ing in the tensor also need to be determined (e.g. the parameter a appearing in 
the elliptic-hyperbolic tensor listed in (j4.10j0 . Once A/, 5 l and all essential pa- 
rameters are known, the transformation from Cartesian coordinates x % to separable 
coordinates u % is 

x* = X/T'(u k ) + 5\ (7.2) 

where x l = T l {u^) is the standard coordinate transformation associated with the 
separable coordinates tabulated in (j4.7j) - (j4.17|) . 

To carry out this procedure, we use the transformation rules relating the param- 
eter matrices of to those of its canonical form (see (J3.6J0 . In matrix form, 
these transformation rules read 

A = X l AX + 2 SiX'Bn) + /x'C^, (7.3a) 
B = X*BX + fj*CX, (7.3b) 
C = A*CA, (7.3c) 

where x = Xx + S, S denotes the symmetric part and 



x = \y \ = x\ x = \y \ = ~x\ X = A x 2 A 2 2 A 3 2 = A 



\7 ) 



[I = | Ai — Ai d" A 2 — A 2 A 3 U <T — A 3 U <T \ = fij 
1x2 \ 2rl \ 1x2 



The identities 

M V = I -5 3 5 1 I , fjtfj* = I 






8 3 


-5 2 \ 




-5 3 









5 2 


-5 1 


o y 





Ax 1 


A 2 X A3 1 




A, 2 


A 2 2 A3 2 


^ ij 


W 


A 2 3 A3 3 






A 2 3 5 2 


A 3 2 5 3 - 


A 3 3 ^ 


V5 3 


A 3 3 5 1 - 


A3 3 * 1 


A 2 2 ^ 




A 3 2 5 1 










) 2 + (5 1 ) 2 






-S 2 S 3 





_ 5 2 5 3 (J 1)2 + (£2)2^ 

(7.4) 

shall prove useful and, in addition, the inverse of ()7.3|) which reads 

A = AAA* + 2 S(XB^) + /xCy , (7.5a) 
£? = ABA* + nCX\ (7.5b) 
C = XCX 1 . (7.5c) 



Our procedure for determining the transformation to canonical form for the 
cases of translational and rotational CKTs is provided in subsections 17.11 and 17.21 
respectively. The three asymmetric CKTs are each treated separately in subsec- 
tions OEHl 

7.1 The translational Killing tensors 

Let us consider first the Cartesian CKT. Since it is necessarily a constant tensor, 
the transformation rules ()7.3j) reduce to A = A*AA, where A = diag(a 1; a 2 , a 3 ), on 
account of (j4.18|) . Trivially, the a« are the eigenvalues of A, the columns of A are 
the (normalized) eigenvectors of A and the translation S is arbitrary. 
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Suppose now that the CKT is circular cylindrical, parabolic cylindrical or elliptic- 
hyperbolic. We may assume without loss of generality that the tensor has the form 
()4.3U|) . since it must necessarily be of this form in order to carry out the classification 
scheme in section |B] Consequently, the rotation and translation are of the form 




— sin0 


°\ 






COS (j) 












U 







and thus the problem reduces to finding <ft, 5 l and 5 2 . The derivation of these 
parameters for each of the three translational CKTs under consideration is provided 
in ^H] (see table 1, p 1432). We now restate these results in our notation. 

1. Circular cylindrical case: The rotation angle <fi is arbitrary and 

£l _ ^23 £ 2 = ^13 

c 3 ' c 3 ' 

2. Parabolic cylindrical case: If 623 7^ 0, then 

bis 

tan</> = — - — , 
023 

and = f for 6 2 3 = (unique mod ir). The components of the translation are 



6 23 (a 2 - ai) + 2a 3 6i 3 , 2 ^13(^2 - ai) - 2a 3 6 



= —r, — 7, ; — ^rr , 
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2(&13 2 + &23 2 ) ' 2(6 13 2 + 6 2 3 2 ) ' 

3. Elliptic-hyperbolic case: Let 

o\ = h 3 2 - 6 2 3 2 + c 3 (a 2 - ai), a 2 = a 3 c 3 - &13&23, A = a 2 + 4a : 
(and note that A is one of the fundamental invariants (|5.7j) ). Then, 

0, if o"2 = and <j\ < 0, 

tan0 = ( 00, if cr 2 = and <j\ > 0, 



2 



£1 

2(T2 



(0 unique mod 7r), and 5 1 and 5 2 are the same as in the circular cylindrical 
case. Moreover, the essential parameter a satisfies 



2 ai — a 2 
a = — 



c 3 c 3 2 



7.2 The rotational Killing tensors 

By the same reasoning used in the previous subsection, we may assume without loss 
of generality that the rotational Killing tensor has the form (j4.31|) . As the isometry 
group for this subspace of rotational webs is the group of translations about the 
z-axis, we set X/ = 5/ and 5 1 = 5 2 = in the transformation rules ([7.3)1 . The 
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determination of 5 3 and hence the transformation to canonical form thus becomes 
a trivial calculation. It follows that 

C 2 

for the spherical, prolate spheroidal and oblate spheroidal cases and 

^3 _ Qi ~ Q3 
~ 2b 12 

for the parabolic case. Finally, the essential parameter a appearing in the transfor- 
mation from Cartesian to prolate (oblate) spheroidal coordinates satisfies 

2 ,a 3 -ai A 2 
a = ± — = = ± XT' 

where Aj and A 2 are the fundamental J(M)-invariants (|5.9|) and the positive (neg- 
ative) signs correspond to the prolate (oblate) spheroidal tensor. 

7.3 The conical case 

The parameter matrices associated with the canonical form of the conical Killing 
tensor specialize to 

A = ail, B = 0, C = diag(ci, c 2 , c 3 ) 

(see ()4.26j) ). From (|7.3cj) . C = X t CX, hence, as in the Cartesian case, the C, are the 
eigenvalues of C and the columns of A are the (normalized) eigenvectors of C. It 
follows that the essential parameters b and c satisfy 

b 2 _ c 2 - ci 

C 2 C 3 - Ci 

thus, in order to satisfy the condition b 2 < c 2 , we can order the eigenvalues such 
that c\ < c 2 < c 3 . Finally, substituting (|7.3cj) into ()7.5b|) leads to B — /xA*C, which 
can easily be solved for d noting the identity (|7.4j) . 

7.4 The paraboloidal case 

The parameter matrices associated with the canonical form of the paraboloidal 
Killing tensor specialize to 

A = diag(Si, a 2 , a 3 ), B = l~b 21 , C = diag(0, 0, c 3 ), (7.6) 

V 0/ 

subject to the constraint 

bwibuhl + c 3 (a 2 - S3)] + hilbuhi + c 3 (Si - S3)] = (7.7) 
(see (14.27)) ). It follows that the essential constants b and c satisfy 



< 



c = 


02—03 

2&12 ' 


c — 


6 = 


(12— Ol 

_26 12 ' 


if 


£3=0, 






b = 


2i— S2 

2bl2 ' 


c — 


b = 


612 
2c 3 ' 


if 


03 7^0, 


621 


= 0, 


c = 


01— S2 
2&21 ' 


c — 


b = 


621 

2 C 3 '_ 


if 


53 7^0, 


bi2 


= 0, 


b = 


2i— S3 
2612 ' 


c — 


b = 


frl2+£>21 
2c 3 ' 


if 


53 7^0, 


b~i2 


^ 0, 62! + 0, 



(7.8) 
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together with the condition b > c. From ()7.3cj) . C = X t CX, hence it follows from 
()7.6|) that C necessarily has a zero eigenvalue of multiplicity two and one other 
eigenvalue c 3 . We now consider the two cases c 3 = and c 3 7^ separately 

If c 3 = 0, then it follows from ()7.7|) that 621 = — &12- Moreover, (j7.3cj) is trivially 
satisfied and (J7.3bj) reduces to B = \ t B\. This implies that B = X l B 2 \, where 
B = diag(— 612 2 , — &12 2 , 0). Therefore, the negative eigenvalue of B 2 determines 
612; we can take its sign to be positive without loss of generality. The normalized 
eigenvectors of B 2 determine A up to a rotation in the eigenspace associated with 
the negative eigenvalue which fixes A up to a parameter if). Finally, it follows that 
(|7.5aj) reduces to A = 2S(BXfi t ) + AAA*. This equation can be solved for the a i: 
6\ and if) which, in general, yields multiple solutions; a particular solution satisfying 
the condition b > c (see (17. 8|) ) can be selected. 

Finally, if c 3 7^ 0, then the eigenproblem ()7.3c|) uniquely determines £3 and A up 
to a parameter ip. Equations (|7.5ajl and (|7.5b|) can then be solved for b\ 2 , b 2 i, 5>i, S\ 
and if), in conjunction with the condition b > c. 

7.5 The ellipsoidal case 

The parameter matrices associated with the canonical form of the ellipsoidal Killing 
tensor specialize to 

A = diag(&i, a 2 , a 3 ), B = 0, C = diag(ci, c 2 , c 3 ), (7.9) 

subject to the constraint 

(Si - a 2 )cid 2 + (o 2 - a 3 )c 2 c 3 + (a 3 - a 1 )c 3 ci = (7.10) 

(see (I4.28jl ). It follows that the essential constants a, b and c satisfy 

a _ b= a^2 y c _ a= 23__2L ; ifg 2 ^ ,£ 3 ^0, 

a-b = ^, c -b=~^, if £2 = 0, c 37 ^0, (7.11) 

together with the condition a > b > c. As in the conical case, (|7.3cJ) implies that the 
5j are the eigenvalues of C and the columns of A are the (normalized) eigenvectors 
of C. There are two cases to consider: (1) the q are all distinct and (2) the q are 
all equal 12 . 

If the Q are all distinct, then the matrix A is uniquely determined (up to the 
ordering of the columns) and 8 can be computed as in the conical case. Upon 
obtaining A and C, the condition a > b > c should be verified and, if necessary, the 
eigenvalues Cj may need to be reordered accordingly. 

If all of the Cj are equal, then ()7.3cj) is trivially satisfied. Equation (J7.5bj) reduces 
to B = ci^tA* which can be solved for <5 using the identity (|7.4jl . Finally, ()7.5aj) 
simplifies to A — ci/x/^* = AAA*. This eigenproblem can be solved to obtain A. 



12 The case of only two equal Ci is impossible, for such a Killing tensor would characterize either 
an elliptic- hyperbolic, prolate spheroidal or oblate spheroidal web. 
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8 Main algorithm 



The above-presented considerations lead to a systematic and computationally effi- 
cient method of determining separable coordinates for the natural Hamiltonian (j4.1|) . 
We emphasize that our algorithm is purely algebraic, and hence is well suited for 
implementation in a symbolic computer algebra system. Indeed, as we mentioned 
in section Q the algorithm has been fully implemented into Maple through the 
KillingTensor package. We now summarize the three main steps of the algorithm. 

(1) Impose the compatibility condition. Using the given potential V in terms of 
Cartesian coordinates x l and a generic Killing tensor K of the form ()3.7jl . impose 
the compatibility condition ([4.3)1 to obtain the equivalent conditions on the Killing 
tensor parameters ([3.9)1 . Computationally, this step amounts to solving a system of 
linear equations in the parameters ([3.9)1 . 

(2) Extract the orthogonal coordinate webs. Decompose the general solution 
obtained in step (1) into the form 

K = e o9 + e 1 K 1 + --- + e n K n , (s.i) 

where li, i = 1, . . . , n are arbitrary constants, g is the metric tensor and {K%, . . . , K n } 
is a linearly independent set of Killing tensors, noting that n ^ 19 since dim /C 2 (E 3 ) = 
20. By theorem 4.1, each Ki must necessarily have normal eigenvectors and distinct 
eigenvalues if it is to characterize separation in one of the eleven separable coordi- 
nate systems. The former can be verified using the TSN conditions ([4.4)1 while the 
latter can be verified efficiently by computing the discriminant of the characteristic 
polynomial of Ki and checking that it does not vanish identically. Finally, relabel 
the K i so that Ki, . . . , K m , m ^ n, are CKTs. 

(3) Classify each Killing tensor and transform to canonical form. For each Ki, 
i = 1, . . . , m, in step (2), classify Ki using the scheme in sectional Finally, using the 
techniques described in section [71 determine the transformation ([7.1)1 which brings 
Ki to its appropriate canonical form. The transformation to separable coordinates 
can be carried out using equation ([7.2)1 . 

Remark. Because of the non-linearity of the integrable eigenvector and distinct 
eigenvalue conditions as well as the non-linearity of the fundamental invariants, cer- 
tain linear combinations of the Ki, i = 1, ... ,n in step (2) may produce Killing 
tensors which characterize separability in coordinate systems not characterized by 
the Ki, i = 1, . . . ,m. In fact, it may be possible to construct such a linear combi- 
nation where the individual Killing tensors of the combination fail to have normal 
eigenvectors or distinct eigenvalues. This will be illustrated by the example in the 
next section. 



9 Application: The Calogero-Moser system 

We apply the algorithm of section |H1 to the three-body inverse square Calogero- 
Moser system with equal masses. It is defined by the natural Hamiltonian ([4.1)1 
with potential 

V = + + . (9.1) 

[x — y) 2 (y — z) 2 (z — x) 2 

Solving the compatibility condition ([4.3)1 with the potential ([9.1)1 yields 

K = a l9 + a 1 K 1 + b 32 K 2 + c 3 K 3 + 73 K 4 , (9.2) 
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where 




— 2yz (x + y — z)z (z + x — y)y^ 
x + y — z)z —2zx (z + y — x)x 
(z + x — y)y (z + y — x)x —2xy 

f y 2 + z 2 —xy —zx 
K$ = ( — xy z 2 + x 2 —yz 
—zx —yz x 2 + y 2 

(9.3) 

Using (|4.4jl . we find that K has normal eigenvectors for all a±, 632, C3 and 73. 
However, it follows that only K 2 and K4 have distinct eigenvalues. It is also useful 
to note that K in (|9.2|) admits a Killing vector 

V = (y- z)X 1 + (z- x)X 2 + (x- y)X 3 . (9.4) 

Using the Killing vector classification scheme in section (jfi.lj) . it follows from ta- 
ble (16. lj) that V is rotational and the transformation 

2 V2> 

./- = A/F +5*, A/ = -)= I -1 y/3 y/2 I , <f = 0, (9.5) 




-1 -v^ v^y 

brings ()9.4|) to the canonical form ()6.3|) . Let us now apply step (3) of the algorithm 
in section |H1 to each of the CKTs in (|9.3J) . 

Applying the transformation (|9.5jl to K 2 yields 

/ 2z -x s 
K% = V3 25 -y 
y— x — y 

It follows that K 2 G /C|.(E 3 ) with respect to the transformed Cartesian coordinates 
x l , and thus, by table EJ we see that K 2 characterizes a parabolic web. Finally, we 
observe from (|4.25|) that K 2 is already in canonical form. Therefore, it follows from 
()9.5|) together with (|4.14j) that the transformation to separable parabolic coordinates 
(/1, u, if)) is given by 

x = ^^ c ° s ^ + - 

y = -^ g /ii/cos^ + ^ 5 /iz/sin^ + ^7 3 (/i 2 - 
z = -^nvcosip - ■j^fiusmip + ^{n 2 

Proceeding similarly for K \. it follows that it enjoys the form 

(2y 2 — z 2 —2xy zx 
-2xy 2x 2 - z 2 yz 
zx yz —x 2 — y 2 

under the transformation ()9.5j) . We conclude from table |U that K % 1 characterizes 
a spherical web and is in canonical form upon comparison with ([4.22)1 . Thus, the 
transformation to separable spherical coordinates is given by ()9.5j) together with 

dnu). 



u 2 ) 



2 2\ 

V ) 
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^5,6 



As mentioned in the remark at the end of section |Hl we can take various linear 
combinations of K±, . . . , in an attempt to find additional separable coordinate 
systems. For example, let K 5 6 = K 3 ±Ki. It follows that K 5 6 both have distinct 
eigenvalues (even though Ki does not). Applying the transformation (|9.5|) to K 5) $ 
yields 

/= Fl + y 2 + z 2 —xy —zx 
—xy =Fl + z 2 + x 2 —yz 
—zx —zy ±2 + x 2 + y 2 

^From table 0] we conclude that K§ and K§ characterize the prolate spheroidal and 
oblate spheroidal webs, respectively. Moreover, K^ 6 is in canonical form (compare 
with (|4.24jl and (|4.25jl ) and the essential parameter appearing in the transformation 
from Cartesian to prolate (or oblate) spheroidal coordinates is a = v3 for both cases 
(see section IT2J) . 

Finally, it follows that Kj = K\ + K 3 + K4 has distinct eigenvalues and admits 
a translational Killing vector 



V = X 1 + x 2 



X 



(9.6) 



It follows that the transformation (|9.5j) brings (|9.6jl to the canonical form (|fj.2j) . 
Applying this transformation to JCf yields 



Ki 



■1 j 



-1 + 3y 2 —3xy 
—3xy — 1 + 3x 2 




which characterizes the circular cylindrical web. 

To summarize, we have shown using the algorithm of section [Sjthat the Calogero- 
Moser system with potential (j9.1|) separates in five orthogonally separable coordinate 
systems, namely, circular cylindrical, spherical, prolate spheroidal, oblate spheroidal 
and parabolic. This result is consistent with that found in [3] and [Oj. 

We conclude by making two remarks. Firstly, our analysis is exhaustive in the 
sense that we have found all possible orthogonally separable coordinate systems 
for which the Calogero-Moser system separates, since it follows that the general 
Killing tensor ()9.2j) admits the rotational Killing vector ()9.4j) and hence can only 
characterize a rotational web. Secondly, the same conclusions from this section hold 
for a weighted Calogero-Moser system with unequal masses. More precisely, the 
Hamiltonian system (|4.1j) with potential 



V 



9i 



+ 



92 



+ 



9?> 



{m\x — m 2 y) 2 {yn2V ~ m 3 z) 2 (m 3 z — m\x) 2 



separates in the five aforementioned coordinate systems, where gt and rrii are con- 
stants and rrii > 0. Moreover, in all five cases, the transformation to separable 
coordinates is given by 



x 



1 - X 1 l T j (u k ) 



+ 8\ \a 



m 1 MN~ 1 
-m 2 m 3 2 MN 
-m 3 m2 2 MN 



m2m 3 M\ 
WL2N m 3 m\M 
—m 3 N m^M/ 



5 l = 0, 



Ki 



where 



M 



f 2 2 1 2 2 1 2 2\-l/2 A7 I 2 1 2\-l/2 

(mi m 2 + m 2 m 3 + m 3 vti\ ) 1 , JM = [1712 + m 3 ) 1 . 



The authors believe that this result is new. 
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10 Conclusions 



In this paper we solve a non-trivial problem of the geometry of orthogonal coordinate 
webs, namely the classification of the eleven orthogonal coordinate webs in E 3 in 
terms of the invariants of the corresponding vector spaces of Killing two-tensors 
and vectors. Notably, the original solution presented here fits well the approach to 
geometry of Felix Klein presented in his Erlangen Program. Moreover, the results 
are successfully applied to the integrability problem of the Hamiltonian systems 
defined in E 3 . From this viewpoint, the well-known Calogero-Moser super-separable 
Hamiltonian system has been integrated within the framework of the (orthogonal) 
Hamilton- Jacobi theory of separation of variables. The other three-dimensional 
pseudo-Riemannian flat space that is amendable to the methods developed in this 
work is Minkowski space E 2 ' 1 . The work in this direction is underway. 
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